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SUMMARY 


Satellite tracking data has been analyzed using three different sets of coef- 
ficients in the mathematical model that describes the earth's gravitational field. 
The results were not intended to be used as a definitive evaluation of the differ- 
ent coefficients but as an assessment of the effects that different sets of coeffi- 
cients, station coordinates, and earth parameters that have been published, have 
on orbital geodetic results. 

Orbital solutions were estimated from optical tracking data from the 
GEOS-I satellite, taken by five major geodetic optical tracking networks. The 
networks and camera types consisted of the SAO Baker-Nunn, GSFC STADAN 
and SPEOPT MOTS 40" and 24", USAF PC-1000, and the US C&GS BC-4. The 
three sets of gravity coefficients used were the SAO M-l set (modified by the 
GEOS-I resonant harmonics), APL 3.5 set, and the NWL 5E-6 set. The semi- 
major axis, gravitational constant, and flattening consistent with each set of co- 
efficients were also used. The station coordinates used were referenced to the 
SAO C-5 standard earth as no other complete set of optical station coordinates 
were available. 

Several long arc orbital analyses were completed using each set of coeffi- 
cients and the results were compared. Orbits were fitted to two overlapping 
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data sets; the arc lengths of these orbits were 5-1/4 days and lday. The orbital 
solutions obtained with each set of coefficients were compared. Furthermore, 
the trajectory differences were computed, and were resolved into radial, cross 
track and along track components . The along track differences were as great 
as 400 meters for the 5-1/4 day arc and 200 meters for the 1 day arc. 

The range measurements of the Goddard Range and Range Rate (GRARR) 
S-Band Tracking System at Rosman were evaluated for 15 passes recorded 
during the first week in January 1966. The actual measurements were com- 
pared with values computed from the optical orbital solutions previously men- 
tioned. The residual differences between the observed and computed ranges 
were analyzed for zero set range bias errors, timing errors, and random 
errors. The error estimates obtained from the orbital solutions determined 
using the SAO M-l gravitational coefficients displayed a much greater degree of 
consistency from pass to pass. Also, the error estimates obtained from the 
5-1/4 day arc were in good agreement with those obtained from the 1 day arc 
when the SAO M-l coefficients were used but agreement was generally poor 
when the reference orbits were determined using the NWL 5E-6 and APL 3.5 
gravitational coefficients. 

Two estimates of the coordinates for the Goddard Range and Range Rate 
station in Tananarive, Madagascar were obtained from independent data sets 
using each set of coefficients. Only the SAO M-l set produced two estimates 
that were consistent; they differed by only 5 meters. 

These comparisons serve to reinforce what was intuitively obvious — that 
for long arc geodetic work, the most complete set of gravity coefficients to- 
gether with consistent station coordinates should be used. 
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GRAVITY MODEL COMPARISON 
USING GEOS-I LONG ARC 
ORBITAL SOLUTIONS 


1.0 INTRODUCTION 


This report presents results and comparisons that have been obtained from 
the reduction of satellite tracking data using three different sets of coefficients 
in the mathematical model that describes the earth’s gravitational field. These 
comparisons were not intended as an evaluation of the coefficients but as an 
assessment of the effects of using the different sets in order to choose the most 
suitable available set of gravity coefficients and station coordinates for long arc 
(greater than 6 revolutions) geodetic purposes. 


These results were obtained using the orbit determination program NONAME 
(Reference 1 ), and the orbital solutions were estimated from optical tracking 
data taken from the GEOS-I satellite. The NONAME program uses a mathe- 
matical function based on Legendre polynomials to approximate the earth's 
gravitational field (Appendix A). Several sets of coefficients for these polynomi- 
als have been published; three of these sets were used for this work; they are: 


1 . The SAO M-l Set. Modified by the GEOS-I Resonant Harmonics (C 13 12 , 
^14,12 * ^15,12 * ® 13, i2’ S i4 1 12 » ^15,12^ (Reference 2 ) , (Reference 7 ) . 

2 . The APL 3.5 Set (Reference 3 ). 


3 . The NWL 5 E -6 Set (Reference 4 ). 


These are presented in Table II. 


For the purposes of these comparisons, the earth's semi -major axis, gra- 
vitational constant, and flattening coefficient that are consistent with each set of 
coefficients were used. These are summarized in Table I. The station coordi- 
nates were unchanged and were referenced to the SAO C -5 standard earth (Ap- 
pendix C), since no other complete set of coordinates for the optical tracking 
stations was available. 


The use of only the one set of station coordinates prevents these results 
being used as any sort of definitive evaluation of these sets of gravity coeffi- 
cients. It should be noted, however, that the ellipsoids defined by the parameters 
in Table I are very similar; thus the station coordinates, if they are fairly 
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Table I 

Parameters for the Earth's Ellipsoid 


Model 

Parameter\^^ 

Gravitational 

Constant 

(km 3 /sec 2 ) 

Semi-Major Axis 
(km) 


Flattening 


SAO M-l 


6378.165 

T. 

298.25 


APL 3.5 


NWL 5E-6 


3.986032 x 10 5 3.986075 x 10 s 3.9860542 x 10 s 


6378.166 

1 . 

298.30 


6378.165 

lT 

298.25 


■ 


Table II 

Harmonic Coefficients (Normalized) 


SAO Ml 


APL 3.5 


2 1 


3 0 

3 1 

3 2 

3 3 


NWL 5E-6 


4 2 


C X 10 6 

S x 10 6 

C x 10 6 

S x 10 6 

C x 10 6 

S x 10 

•484.1735 


-484.198 


-484.194 






0.016 

0.062 

2.379 

-1.351 

2.381 

-1.198 

2.446 

-1.511 

0.9623 


1.011 


0.984 


1.936 

0.266 

1.84 

0.215 

2.148 

0.27^ 

0.734 

-0.538 

1.219 

-0.6791 

0.978 

-0.90( 

0.561 

1.620 

0.6609 

0.9795 

0.585 

1.62E 

0.5497 


0.467 


0.507 


-0.572 

-0.469 

-0.5624 

-0.4403 

-0.495 

—0.57? 

0.330 

0.661 

-0.4179 

0.4438 

0.274 

0.671 




Table II (Continued) 



SAO Ml 

APL 3.5 

NWL 5E-6 

C x 10 6 S x 10 6 

C X 10 6 S X 10 6 

C x 10 6 S x 10 6 


0.851 -0.190 


-0.053 


0.230 


0.8464 

-0.2106 


0.007062 

0.1898 


1.030 —0.247 

-0.413 0.33( 


•0.0633 

-0.079 

0.631 

-0.520 

•0.265 

0.156 


-0.103 

A o on 

0.007 

0.064 

-0.592 


0.084 

0.1370 

0.2684 

0.09131 

-0.4884 

-0.03358 


-0.1669 

-0.3379 

0.1035 

-0.260 

- 0.6686 


0.045 

0.032 

0.637 

-0.389 

-0.549 

0.215 


-0.11J 

-0.325 

-0.124 

0.145 

-0.594 


-0.1792 

-0.047 

0.069 

-0.054 

-0.044 

-0.313 

-0.040 


-0.027 

-0.366 

0.031 

-0.518 

-0.458 

-0.155 


-0.103 

-0.0002093 

-0.1610 

0.5303 

-0.3069 

-0.18 

0.01434 


0.1009 

-0.1555 

0.05111 

-0.5087 

-0.5091 

-0.2316 


-0.219 

-0.085 

0.129 

- 0.020 

-0.193 

-0.093 

-0.324 


0.19S 

-0.451 

-0.134 

-0.31( 

-0.78( 

-0136 


-0.0860 

0.197 

0.364 

0.250 

-0.152 

0.076 

-0.328 



0.153 

0.156 

0.1261 

0.163 

0.4586 

0.018 

0.3938 

-0.102 

-0.1368 

0.054 

-0.05682 

0.063 

-0.4552 


0.09355 

0.05998 

-0.2067 

0.0004798 

-0.1871 

0.758 


0.105 

0.331 

0.350 

0.323 

-0.467 

0.055 

-0.477 


0.055 


0.096 


0.08840 


-0.1443 




Table II (Continued) 


n 


SAO Ml 

APL 3.5 

NWL 5E-6 

m 

C x 10 6 

S x 10 6 

ni 

X 

H- 1 

O 

S x 10 6 

C x 10 6 

S x 10 6 

8 

0 

0.0655 


0.170 




8 

1 

-0.075 

0.065 

-0.1481 

-0.04843 



8 

2 

0.026 

0.039 

0.09472 

-0.03764 



8 

3 

-0.037 

0.004 

-0.05497 

0.2168 



8 

4 

-0.212 

-0.012 

-0.06901 

0.03761 



8 

5 

-0.053 

0.118 

0.08040 

-0.002495 



8 

6 

-0.017 

0.318 

-0.02193 

0.6658 



8 

7 

-0.0087 

0.031 

0.1697 

-0.07009 



8 

8 

-0.248 

0.102 

-0.1457 

0.09424 



9 

0 

0.0122 


0.041 




9 

1 

0.117 

0.012 





10 

00 

0.0118 






10 

01 

0.105 

-0.126 





10 

02 

-0.105 

-0.042 





10 

03 

-0.065 

0.030 





10 

04 

-0.074 

-0.111 





11 

00 

-0.0630 


0.104 




11 

01 

-0.053 

0.015 





12 

00 

0.0714 


0.062 




12 

01 

-0.163 

-0.071 





12 

02 

-0.103 

-0.0051 





12 

12 

-0.031 

0.0008 
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Table II (Continued) 


n 

m 

SAO Ml 

APL 3.5 

NWL 5E-6 

C x 10 6 

S x 10 6 

C X 10 6 

S x 10 6 

C x 10 6 

S x 10 6 

13 

00 

0.0219 






13 

12 

-0.06769 

0.06245 





13 

13 

-0.059 

0.077 

-0.4689 

0.04748 

-0.03 

0.11 

14 

00 

-0.0332 






14 

01 

-0.015 

0.0053 





14 

11 

0.0002 

-0.0001 





14 

12 

0.00261 

-0.02457 





14 

14 

-0.014 

-0.003 

-0.06368 

-0.037 



15 

09 

-0.0009 

-0.0018 





15 

12 

-0.07473 

-0.01026 





15 

13 

-0.058 

-0.046 



-0.06 

-0.06 

15 

14 

0.0043 

-0.0211 

0.00087843 

-0.0101 

0.01 

-0.03 


accurately determined with reference to the center of mass, as the SAO C-5 
coordinates are generally accepted to be, should not introduce any large differ- 
ences in the results. 

Several long arc analyses were completed using each set of coefficients in 
turn, and the results have been compared; these are discussed in some detail in 
Sections 3.0 - 6.0. 


2.0 DESCRIPTION OF THE EARTH'S GRAVITATIONAL FIELD 

The earth's geopotential can be approximated by the following mathe- 
matical model: 


5 
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GM 


( 1 ) 


f k n 



(sin 


*)[ 


C cos mk + S sin m k 




where 

G is the universal gravitational constant, 

M is the mass of the earth, 
r is the geocentric satellite distance, 
a is the earth's mean equatorial radius, 

4 > is the sub-satellite geocentric latitude, 

A. is the sub-satellite east longitude, 

p™ (sin 4 >) are the associated Legendre polynomials of degree n and order m, 

and 

C nm , S nm are the denormal ized gravitational coefficients. 

The denormalized gravitational coefficients are related to the normalized 
coefficients (C , S ) as indicated below: 

C(n,m) = ^(n ~ m) ! (2n + 1) K/ (n + m) ! J A C(n,m) 

S(n,m) = [(n - m) ! (2n + 1) K/(n + m) ! ] A S(n,m) 

where , 

K = 1 when m = 0 
K = 2 when m ^ 0 

The geopotential formulated in this manner can be converted into gravita- 
tional accelerations in inertial coordinates (x, y, z) as follows: 
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where the subscript "o" denotes accelerations due to the earth's gravitational 
field. Similar expressions hold for y m and z e . The NONAME program uses a 
model in this form to compute the accelerations due to the earth's gravitational 
field. 

The three different sets of harmonic coefficients (normalized) and associated 
earth parameters used in this analysis are shown in Tables I and II. The SAO M-l 
is the largest set with a total of 122 coefficients; the APL 3.5 set has 84 coeffi- 
cients and the NWL 5E-6 set has 64. Of these three sets, only the SAO M-l set 
has GEOS-I resonant terms (harmonics of order 12). 

The SAO M-l set was determined using optical observations from a number 
of satellites, and the other two sets were determined from Tranet Doppler obser- 
vations, again from a number of satellites. 


3.0 DIFFERENCES IN ORBITAL SOLUTIONS 

Orbits were fitted to two data sets from the first week in January 1966, and 
the arc lengths of these orbits were 5-1/4 days and 1 day. The 5-1/4 day arc 
covered the period from 01 hrs. GMT on December 31, 1965 to 06 hrs., Janu- 
ary 5, 1966, and the data set consisted of 1057 optical observations.* The 1 day 
arc covered the period from 06 hrs., January 2, 1966, to 08 hrs., January 3, 1966. 
This data was a subset of the 5-1/4 day arc data set and consisted of 444 optical 
observations. These data sets are summarized in Table III. 

The root mean squares of the observations about the orbital solutions were 
computed and these are shown in Table IV. The r.m.s. values were lower for the 
orbits fitted using the SAO M-l set for both arcs. The differences between the 
observed measurements and values computed from the orbital solutions were 
computed and plotted on histograms; these are shown in Figures 1-4. The right 
ascension residuals shown in the Figures have been multiplied by the cosine of 
the corresponding declination measurement to account for the degradation of the 
measurements recorded when the declination value was large. These Figures 
clearly indicate that the orbital solutions obtained with the SAO M-l set of coeffi- 
cients fit the data sets better than the other solutions. This is especially true 


* 


Right ascension plus declination measurements. 
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Table in 

Summary of Optical Measurements By Station 


Network 

Station 

Camera Type 

No. of Observations 

5- 1/ 4 Day Arc 

1 Day Arc 

SAO 

10RGAN 

Baker -Nunn 

2 



1JUPTR 

IT 

26 

26 


1NATOL 

IT 

8 

2 


OSLONR 

H 

4 



AUSBAK 


4 



1SHRAZ 

TT 

2 

2 


1SPAIN 

TT 

6 



1 TOKYO 

M 

12 

4 


1VILDO 

ff 

2 



1MAUIO 

IT 

2 



AGASSI 

Geodetic 36" 

10 



Total: 


78 

34 

SPEOPT 

1COLBA 

MOTS 40" 

164 

71 


1JUM40 

TT 

22 

16 


1BERMD 

TT 

84 

36 


1PURI0 

II 

14 



1DENVR 

T! 

70 

14 


1JUM24 

MOTS 24" 

26 

24 


Total: 


380 

158 

STADAN 

1FTMYR 

MOTS 40" 

82 

54 


1BPOIN 

M 

53 



1GFORK 

TT 

26 

9 


1MOJAV 

II 

25 

25 


Total: 


186 

91 

USAF 

HUNTER 

PC-1000 

59 

47 


SWANIS 

TT 

14 

14 


GRDTRK 

IT 

7 



ANTIGA 

M 

26 



SEMMES 

TT 

60 

36 


CURACO 

TT 

40 

26 


HOME ST 

TT 

94 

24 


JUPRAF 

TT 

17 

17 


BEDFRD 

TT 

22 



ABERDN 

TT 

74 



Total: 


413 

164 

Total of All Observations 

1,057 

444 
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- 10 -8 -6 -4 -2 0 2 4 6 
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SECONDS OF ARC 

Figure 2. Declination Residuals from 1 Day Arc 





FREQUENCY FREQUENCY FREQUENCY 


SAO Ml GRAVITY 



SECONDS OF ARC 

Figure 4. Declination Residuals from 5 1/4 Day Arc 
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Table IV 

Root Mean Squares About The Orbital Solution 


Arc Length 

R.M.S. of Fit (Secs, of Arc) 

SAO M-l 

APL 3.5 

NWL 5E-6 

5-1/4 day 

3.08 

11.14 

11.01 

1 day 

2.33 

2.50 

4.54 


for the 5-1/4 day arc (Figures 3-4) where the distributions of the right ascension 
and declination residuals for the APL 3,5 and NWL 5E-6 solutions are obviously 
non normal. 


4.0 TRAJECTORY DIFFERENCES 

Satellite position differences were computed at five minute intervals using 
the orbital solution obtained from the SAO M-l coefficients as a standard for 
comparing the solutions obtained for the same data set when the APL 3.5 and 
NWL 5E-6 gravitational coefficients were used. These position differences were 
resolved into along track, cross track and radial components and are shown in 
Figures 5-8. 

The along track differences were the largest. They were as large as 400 
meters for the 5-1/4 day arc and 200 meters for the 1 day arc. The cross track 
and radial differences were approximately the same order of magnitude and 
were as large as 100 meters for the 5-1/4 day arc and 50 meters for the 1 day 
arc. The differences have a period approximately equal to the period of the 
satellite (2 hrs), and, in addition, the along track differences have some other 
long period associated with them. The periods of the along track, cross track, 
and radial differences are not in phase, and in general, the minima occur where 
there is good data coverage; this is shown by the solid blocks in the figures. 


5.0 EVALUATION OF THE ROSMAN GRARR RANGE ACCURACY 


The range measurements of the Goddard Range and Range Rate (GRARR) 

S Band Tracking System at Rosman, North Carolina, were evaluated by comparing 
the actual measurements with values computed from the optical reference orbits 
(Reference 5). The 5-1/4 day arc and the 1 day arc discussed in Section 3.0 were 
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Figure 5. Differences Between Trajectories Obtained from SAO M-l and APL 3.5 Gravity Models 1 Day Arc 
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used as the reference orbits . A summary of the station pass times , simultaneous 
optical coverage and maximum elevation angles for 15 GRARR passes recorded 
at Rosman during the 5-1/4 day optical arc is presented in Table V. It should be 
noted that a large percentage of the optical data used in the determination of the 
reference orbital solutions was recorded by tracking stations located on or near 


Table V 

Summary of GRARR Passes At Rosman Occurring 
During The 5-1/4 Day Optical Orbital Arc 


Pass 

No. 

Transponder 

Channel* 

Date 

Time 

No. of Obs. in Pass 

Max. 

Elevation 

Angle 

R/R 

Optical 

652 

A 

12/31/65 

06 h 

18 

18 

31.3° 

653 

A 

12/31/65 

O 

00 

x 

28 

30 

65.4° 

664 

A 

1/1/66 

06 h 

28 

78 

36.6° 

665 

A 

1/1/66 

08 h 

32 

95 

51.8° 

673 

A 

1/1/66 

23 h 

34 

0 

53.5° 

676 

A 

1/2/66 

06 h 

32 

106 

43.3° 

677 

A 

1/2/66 

08 h 

28 

138 

40.2° 

685 

A 

1/2/66 

23 h 

34 

10 

46.5° 

688 

A 

1/3/66 

06 h 

30 

101 

52.2° 

689 

A 

1/3/66 

X 

00 

o 

14 

79 

30.1° 

697 

A 

1/3/66 

23 h 

44 

0 

40.8° 

700 

C 

1/4/66 

06 h 

36 

100 

62.7° 

708 

C 

1/4/66 

21 h 

48 

0 

84.2° 

709 

C 

1/4/66 

23 h 

42 

14 

35.8° 

712 

A 

1/5/66 

06 h 

36 

66 

76.6° 


*The GEOS-I GRARR transponder contained two channels denoted A and C which received signals 
at 2271.9328 MHz and 2270.1328 MHz respectively. 
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the North American Continent (Table III) and also most of the Rosman GRARR 
passes had simultaneous optical data coverage. 

For each GRARR pass over Rosman, zero set range bias errors, timing 
errors, and random errors were estimated from the residual differences be- 
tween the observed and calculated ranges; these are summarized in Tables VI- 
VIII. The 12 A channel passes had a mean zero set error of -10 meters with a 
standard deviation of 8.8 meters, and a timing error of -2.4 milliseconds with a 
standard deviation of 2.4 milliseconds when compared with the orbital solution 
determined with the SAO M-l gravitational coefficients; a mean zero set error 
of -33.7 meters with a standard deviation of 21.0 meters, and a timing error of 
-0.8 milliseconds with a standard deviation of 9.8 milliseconds when compared 
with the orbital solution determined with the NWL 5E-6 coefficients; and a mean 
zero set error of -34,6 meters with a standard deviation of 30.7 meters, and a 
timing error of 3.1 milliseconds with a standard deviation of 22.1 milliseconds 
when compared with orbital solution determined with the APL 3.5 coefficients. 
As indicated by the standard deviations associated with these errors, the esti- 
mates obtained from the orbital solutions fitted using the SAO M-l set of coeffi- 
cients were significantly less variable than those obtained using the other two 
sets. In addition, the estimates obtained from the shorter overlapping 1 day arc 
were only consistent with the 5-1/4 day arc estimates when the orbital solutions 
were obtained with the SAO M-l set of coefficients. 


6.0 ESTIMATION OF COORDINATES FOR THE GRARR MADGAR SITE 

Two independent estimates of the coordinates of the GRARR site in 
Tananarive, Madagascar (MADGAR) were obtained using each set of coefficients 
(Reference 6). One estimate was obtained from optical flash sequence data re- 
corded by the MOTS 40" camera (1 TANAN) during July 1966 and the other from 
range measurements taken at MADGAR during November 1965. The data sets 
used for these estimations are shown in Tables IX and X. 

The two sets of coordinates estimated using the SAO M-l coefficients were 
very consistent, within 5 meters of each other; whereas the estimates obtained 
using the other two sets of coefficients were not at all consistent. These esti- 
mates are shown in Table XI and Figure 9. 
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Table VI 

Summary of Rosman Zero-Set Range Bias Error Estimates 

(meters) 


Pass No. 


Transponder 

Channel 


SAO Ml NWL 5E-6 API 


No. 1* No. 2 No. 1 No. 2 No. 1 No. 2 



Mean 


Std. dev. 



-51.8 - 23.6 -20 

-60.0 -112.5 -46 



No. 1 and No. 2 refer to the 5 l A day and 1 day orbital arcs described in Section 3.0 














Table Vin 

Summary of Random Error Estimates 
(meters) 


Pass No. 

Transponder 

SAO 

Ml 

NWL 

5E-6 

APL 

3.5 

Channel 

No. 1 



No. 2 


No. 2 

652 

A 

3.3 


3.5 


6.0 


653 

A 

3.6 


10.5 


7.3 


664 

A 

3.2 


3.6 


6.1 


665 

A 

4.5 




3.6 


673 

A 

2.6 


4.7 


3.1 


676 

A 

2.9 

3.0 

3.9 

5.3 


3.8 

677 

A 

3.2 

4.0 

10.4 

4.7 

2.8 

5.3 

685 

A 

2.5 

2.3 

3.5 


2.6 

5.2 

688 

A 

4.5 


8.3 

10.3 

7.0 


689 

A 

4.1 

3.8 

7.4 

5.0 

3.9 

3.9 

697 

A 

3.4 


3.5 


3.8 


700 

C 

5.8 


10.8 


16.8 


708 

C 

4.5 


3.9 


5.3 


709 

c 

3.1 


4.4 


3.6 


712 

A 

6.0 






Mean 

A 

3.7 


6.6 


6.2 


Mean 

C 

4.5 


6.4 


8.6 
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Table IX 


Summary of Optical Data by Station 
For July 9, 10, and 11, 1966 


Station 

No. of Measurements 

Right Ascension 

Declination 

1TANAN 

14 

14 

1ROSMA 

7 

7 

1COLBA 

14 

1 A 

j-r 

1BPOIN 

14 

14 

1DENVR 

20 

20 

1JOBUR 

14 

14 

10RGAN 

91 

91 

lOLFAN 

28 

28 

1SPAIN 

21 

21 

1QUIPA 

28 

28 

1CURAC 

28 

28 

1JUPTR 

35 

35 

1VILDO 

7 

7 

AUSBAK 

14 

14 

1MAUIO 

28 

28 

EDWAFB 

2 

2 

Total 

365 

365 
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Table X 


Summary of Optical and GRARR Data By Station 
For November 28 and 29, 1965 


Station 

No. of Measurements 

Right Ascension 

Declination 

10RGAN 

59 

59 

10LFAN 

1 

1 

ISP AIN 

1 

1 

1QUIPA 

2 

2 

1CURAC 

96 

96 

1JUPTR 

127 

127 

1VILD0 

1 

1 

Total 

287 

287 


Range 


MADGAR 

24 
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Table XI 

Estimated Coordinates for Madgar 


SAP Ml Gravity Model 



Latitude 

E. Longitude 

Spheroid Height 

Optical Estimate 

-19° V 19.5” 

47° 18’ 7.9" 

1380.0 meters 

GRARR Estimate 

-19° 1' 19.4” 

© 

00 

00 

rH 

o 

1382.6 meters 

Difference 

0.1” 

-0.1" 

-2.6 meters 


APL 3.5 Gravity Model 



Latitude 

E. Longitude 

• 

Spheroid Height 

Optical Estimate 

-19° 1' 22.6" 

47° 18’ 5.1" 

1454.5 meters 

GRARR Estimate 

-19° 1' 23.7" 

47° 18’ 5.7" 

1443.1 meters 

Difference 

-1.1" 

-0.6" 

11.4 meters 


NWL 5E-6 Gravity Model 




Latitude 

E. Longitude 

Spheroid Height 

Optical Estimate 

-19° 1' 22.9" 

47° 18’ 

4.8" 

1458.0 meters 

GRARR Estimate 

-19° 1’ 24.9" 

00 

1—1 

o 

C> 

7.5" 

1467.4 meters 

Difference 

-2.0" 


-2.7" 

10.6 meters 
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LATITUDE 


EAST LONGITUDE 


8’ 3.0’ 


4.0” 5.0” 


7.0” 8.0’ 



-19°1’ 18.0’ 


1 OPTICAL ESTIMATE - SAO Ml GRAVITY 

2 GRARR ESTIMATE - SAO Ml GRAVITY 

3 OPTICAL ESTIMATE - APL 3.5 GRAVITY 

4 GRARR ESTIMATE - APL 3.5 GRAVITY 

5 OPTICAL ESTIMATE - NWL 5E-6 GRAVITY 


6 GRARR ESTIMATE 


NWL 5E-6 GRAVITY 


Figure 9. Estimated Coordinates for MADGAR 
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APPENDIX A 


FORCE MODELS USED IN THE 
NONAME ORBIT DETERMINATION SYSTEM 

1.1 FORCE MODELS 

The data reduction program in its present form incorporates four force 
models. These are: 

1. The earth gravitational field 

2. The solar and lunar gravitational perturbations 

3. Solar radiation pressure 

4. Atmospheric drag 

The program is designed such that the gravitational coefficients and pertinent 
physical characteristics of satellites, such as reflectivity, cross-sectional area 
mass, and drag coefficient can be simply changed through card input or block 
data statement. 


1.2 THE EARTH'S GRAVITATIONAL FIELD 


The formulation of the geopotential used is: 


\ Tl (?) P ” (sin 0) [ Cnm C ° S + Snm sinm \| \ (A_1) 


u - 


n= 2 m = 0 


where 


G 

M 


is the universal gravitational constant 

is the mass of the earth 

is the geocentric satellite distance 


A-l 


a 


is the earth's mean equatorial radius 


4> is the sub-satellite geocentric latitude 

k is the sub-satellite east longitude 

p” (sin <fi) is the associated spherical harmonic of degree n and order m. 

The design of the potential function requires that denormalized gravitational co- 
efficients C n m and S n m be used. The program is presently capable of accepting 
coefficients up to (20,20) or any sub-set of these. 

The transformation of the geopotential in earth-fixed coordinates (r, <fi, k) 
to gravitational accelerations in inertial coordinates (x, y, z) is accomplished 
as follows: 

.. _ Bu Br Bu B0 Bu BA. 

Xl ® “ b7b^ + B0B 7 + bTb 7 ; y ® ’ Zl ® (A_2) 

where the subscript "®" denotes accelerations due to the earth's field. 


1.3 SOLAR AND LUNAR GRAVITATIONAL PERTURBATIONS 

The perturbations caused by a third body, e.g., the sun or moon, on a satel- 
lite orbit are treated by defining a disturbing function R (Reference 1) which can 
be treated as the potential function U. For the solar perturbation R 0 takes the 
form 


R 


© 



(A -3) 


where 

S = cos (r, r G ) 

m 0 is the mass of the sun in earth masses 
r 0 is the geocentric position vector of the sun 
r 0 is the geocentric distance to the sun 


A -2 


r is the geocentric distance to the satellite 

r is the geocentric position vector of the satellite 

G is the universal gravitational constant 
M is the mass of the earth 

The acceleration of the satellite due to the sun is then 

_ B r 30 ^R© BA. 

*® Br Bx B0 Bx BA. Bx’ ^® ’ z ® ^ 

where <fi and A. are the latitude and longitude of the satellite respectively. The 
lunar perturbations are found from Equation (A -3) by substituting the lunar mass 
and distance fox' those of the sun. 

The lunar and solar ephemerides are computed internal to the program. 
These positions are computed at ten equal intervals over each five day period 
and least squares fit to a fourth order polynomial in time about the midpoint of 
the five day period. The positions of these bodies are then determined at each 
data point by evaluating the polynomial at the observation time. 


1.4 SOLAR RADIATION PRESSURE 


The acceleration acting on a satellite due to solar radiation pressure is 
formulated as follows (Reference 2). 


AP ° 

X RAD “ “ ~ L , ; >W Z RAD ,A " 5 ' 


where 

L is the inertial unit vector from the geocenter to the sun and whose com- 
ponents are L x , L y) L z 

A is the cross sectional area of the satellite 

m is the satellite mass 

y is a factor depending on the reflective characteristics of the satellite. 
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v is the eclipse factor such that: 


{ 0 when satellite is in earth's shadow 
1 when satellite is illuminated by the sun 

P e is the solar radiation pressure in the vicinity of the earth, 


4.5 x 


10 “ 


Newton 


m 


At present, it is assumed that the satellite is specularly reflecting with re- 
flectivity, p, and thus 


7 - (1 + p) . (A-6) 

The vector L and the eclipse factor are determined from the solar ephemeris 
subroutine previously described, the satellite ephemeris, and involve the approxi- 
mation of a cylindrical earth shadow. 


1.5 ATMOSPHERIC DRAG 

The atmospheric decelerations are computed as follows: 


pC D Avv x 

X DRAG ~ 2ni ’ y DRAG’ Z DRAG 


(A-7) 


where 

p is the ambient atmospheric density 
C D is the satellite drag coefficient 

A is the projected area of the satellite on a plane perpendicular to direction 
of motion 

m is the satellite mass. 

The velocity vector v given in inertial coordinates by 


A-4 


V 


(A -8) 


" + V yi + V z k 

can be chosen to be either the velocity relative to the atmosphere which implies 
that the atmosphere rotates with the earth or the inertial velocity which assumes 
that the atmosphere is static. Presently, the former assumption is made. 

The density, p, is computed from the 1962 U.S. Standard Atmosphere. 
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PREPROCESSING OF OPTICAL OBSERVATIONS 


2.1 PREPROCESSING OF OPTICAL DATA 

The first step in the processing of optical observations is usually performed 
by the observing source. This consists of developing a plate or film, identifying 
the image or images of the satellite and the images of several reference stars 
whose right ascensions and declinations are well known. The initial measure- 
ments of both satellite images and reference stars consist of linear rectangular 
coordinates. From the knowledge of the spherical coordinates of the reference 
stars, the right ascensions and declinations of the satellite images may be calcu- 
lated, These coordinates as received by the preprocessor may be referred to the 
mean equator and equinox of date, true equator and equinox of date, or mean 
equator and equinox of some standard epoch. 

Preprocessing includes, for example in the case of the GEOS-I SAO Baker - 
Nunn data, updating the observations from a mean equinox and equator of 1950.0 
to the true equinox and equator of date through the luni-solar precession and nu- 
tation effects, the correction for planetary aberration, and the transformation of 
the A-S (SAO atomic time) time tag to UTC time. It is necessary to know UT1 
time when the angle between Aries and the Mean Greenwich Meridian is re- 
quired, UT1 time is then calculated on the basis of the differences (UT1-UTC) 
as published by the U.S. Naval Observatory. In the case of active flash data, 
where the time is recoverable to better than 100 microseconds through the use 
of APL published corrections to the satellite on-board clock (Reference 1), the 
time tag is shifted to correspond to the center of the photographic flashing light 
image. This latter adjustment corresponds to a shift of 0.5 milliseconds which 
is equivalent to approximately 4.0 meters of satellite position. 

Currently, the preprocessor transforms all right ascensions and declinations 
to the true equator and equinox of the epoch of the observations being processed. 
If the observations were originally referred to the mean equator and equinox of 
a particular epoch, it is only necessary to precess from that epoch to the dates 
of the observation. However, if they were referred to the true equator and 
equinox of a particular epoch, it is necessary first to transform them to the 
mean equator and equinox of that same epoch and then precess to the epochs of 
the observations. 

Finally, a transformation must be made from the mean equator and equinox 
of the epoch of the observation to the true equator and equinox of the epoch of the 
observation. 
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2.2 NUTATION 


The transformations from the true equator and equinox of date to the mean 
equator and equinox of date is 


Y = NX 


(B-l) 


where 


Y 


cos S m cos a 

m m 

cos S sin a 

m m 


sin S 

m 


x 


cos 8 t cos a T 
cos § T sin a T 

sin S T 


(B-2) 


(B-3) 


1 


+ All/ COS € 
r m 


+ A \p sin e m 


N 


- Ewp cos e m 


1 


+ Ae 


(B-4) 


- &\p sin e m 


- Ae 


1 


where 

a , 8 = right ascension and declination referred to mean equator and 

equinox of date 

a T , 8 X = right ascension and declination referred to true equator and 
equinox of date 


B-2 


e ra = mean obliquity of date 

Ai/- = nutation in longitude 

Ae = nutation in obliquity. 

The inverse transformation is simply: 

X = N" 1 Y = N t Y (B-5) 


2.3 PRECESSION 

The transformation from the mean equator and equinox of 1950.0 to the mean 
equator and equinox of an arbitrary epoch t 1 is 

Y = PX (B-6) 

where 


Y 


cos S tl cos a t x 
cos 5 u sin a (1 
sin 8 x 


X 


COS 5 19 50 , 0 
COS ^1950.0 
S * n ^1950. 0 


cos a 
sin a 


1950. o 
1950.0 


(B-7) 


(B-8) 


B-3 




Since the expression for z, 9, £ are tied to 1950.0 as an epoch, the preces- 
sion between 2 different epochs, neither of which is 1950.0, must be performed 
in two steps, using 1950.0 as an intermediary epoch. 
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TRACKING STATION COORDINATES 


1.0 DATUM PARAMETERS AND STATION COORDINATES 

For the purpose of long arc satellite data reduction and intercomparison, 
all GEOS-I participating tracking stations have been transformed to a common 
datum. The common datum selected is the SAO Standard Earth C-5 Model 
(Reference 1) in which the Baker-Nunn station positions are used as the con- 
trolling stations for all other stations to be transformed. The semi-major axis 
and flattening coefficient for the SAO C-5 Earth Model are 6378165 meters and 
298.25 respectively. Descriptions and formulations to effect the transformations 
from major and isolated datums are presented in Reference 2. The transforma- 
tion of local datum station coordinates to a common center of mass reference 
system is important to be performed since the datum shifts are quite large. 

For example, on the North American datum the center of mass shift to the C-5 
Standard Earth is approximately 250 meters. The center of mass coordinates of 
the SAO C-5 Baker-Nunn stations are assessed by SAO to have approximately 
20 meter accuracy. 

In order to effect any transformation, the parameters of the original datums 
must be known as well as the geodetic latitude, longitude and height. Table I 
provides a listing of the original datums and their parameters on which the sta- 
tions were originally surveyed. Tables II to XI list alternately the original sur- 
veyed ellipsoidal positions and the SAO C-5 ellipsoidal positions for over 100 
GEOS-I tracking stations that have been used in the long arc intercomparison 
effort. These tables contain symbols designating the source of original station 
coordinates. The symbols are defined in Section 2.0 with a list of source infor- 
mation. The C-5 positions for 1TANAN and MADGAR (Reference 3) have been 
derived by the station estimation technique contained in the Orbit Determination 
Program NONAME. Tables XII to XXI provide a listing of the proper station 
names from which the six letter designations have been derived. 
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Table I 


Parameters of Original Datums 


Datum Name 

Semi-Major Axis 
(meters) 

1/f 

North American (N.A.) 

6378206.4 

294.9787 

European 

6378388.0 

297.0 

Tokyo 

6377397.2 

299.1528 

Argentina 

6378388.0 

297.0 

Mercury 

6378166.0 

298.3 

Madagascar 

6378388.0 

297.0 

Australian Nat'l. 

6378160.0 

298.25 

Old Hawaiian 

6378206.4 

294.9787 

Indian 

6377276.3 

300.8017 

Arc (Cape) 

6378249.1 

293.4663 

1966 Canton Astro 

6378388.0 

297.0 

Johnston Island 



1961 

6378388.0 

297.0 

Midway Astro 1961 

6378388.0 

297.0 

Navy Iben Astro 



1947 

6378206.4 

294.9787 

Provisional DOS 

6378388.0 

297.0 

Astro 1962, 65 



Allen Sodano Lt. 

6378388.0 

297.0 

1966 SECOR ASTRO 

6378388.0 

297.0 

Viti Levu 1916 

6378249.1 

293.4663 

CORREGO ALEGRE 

6378206.4 

294.9787 

USGS 1962 ASTRO 

6378206.4 

294.9787 


BERNE 


6377397.2 


299.1528 







Table II 

SAO — Optical — Source A 


— 

Source 

Name 

Station 

No. 

— 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 


10RGAN 

9001 

32°25'24'.'56 

253°26'51"17 

1649 

N.A. 




32 25 24.70 

253 26 48.29 

1610 

C-5 


lOLFAN 

9002 

-25 57 33.85 

28 14 53.91 

1562 

Arc (Cape) 




-25 57 37.67 

28 14 51.45 

1560 

C-5 


WOOMER 

9003 

-31 06 07.26 

136 46 58.70 

185 

Australian 




-31 06 04.14 

136 47 01.93 

158 

C-5 


1SPAIN 

9004 

36 27 51.24 

353 47 41.47 

7 

European 




36 27 46.68 

353 47 36.55 

56 

C-5 


1 TOKYO 

9005 

35 40 11.08 

139 32 28.22 

58 

Tokyo 




35 40 23.03 

139 32 16.42 

84 

C-5 


1NATOL 

9006 

29 21 38.90 

79 27 25.61 

1847 

European 




29 21 34.38 

79 27 27.05 

1855 

C-5 


1 QUIP A ; 

9007 

-16 28 05.09 

288 30 22.84 

2600 

N.A. 


j 


-16 27 58.04 

288 30 24.02 

2479 

C-5 


1SHRAZ 

9008 

29 38 17.96 

52 31 11.80 

1578 

European 


i 


29 38 13.59 

52 31 11.20 

1561 

C-5 


1C UR AC 

9009 

12 05 21.55 

291 09 42.55 

23 

N.A. 




12 05 24.93 

291 09 43.97 

-33 

C-5 


1JUPTR 

9010 

27 01 13.00 

279 53 12.92 

26 

N.A. 




27 01 14.23 

279 53 12.95 

-36 

C-5 


1VILDO 

9011 

-31 56 36.53 

294 53 39.82 

598 

Argentinean 




-31 56 36.35 

294 53 36.11 

636 

C-5 


1MAUIO | 

i 

9012 

20 42 37.49 

203 44 24.11 

3027 

Old Hawaiian 




20 42 25.66 

203 44 33.23 

3027 | 

C-5 


C-3 














Table II (Continued) 


Source 

Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

A 

AUSBAK 

9023 

-31°23'3Cf.'82 

136°52'39’.'02 

164 

Australian 




-31 23 27.69 

136 52 42.23 

137 

C-5 

A 

OSLONR 

9426 

60 12 40.38 

10 45 08.74 

585 

European 




60 12 38.88 

10 45 02.26 

573 

C-5 

I 

NATALB* 

9029 

-05 55 50.00 

324 50 18.00 

112 

N.A. 




-05 55 43.49 

324 50 21.30 

45 

C-5 

D 

AGASSI* 

9050 

42 30 20.97 

288 26 28.71 

193 

N.A. 




42 30 20.51 

288 26 29.79 

138 

C-5 

I 

COLDLK* 

9424 

54 44 38.02 

249 57 25.85 

597 

N.A. 




54 44 37.26 

249 57 21.90 

548 

C-5 

I 

EDWAFB* 

9425 

34 57 50.68 

242 05 11.39 

784 

N.A. 




34 57 50.17 

242 05 07.80 

754 

C-5 

I 

RIGLAT* 

9428 

56 56 54.00 

24 03 42.00 

5 

European 




56 56 52.37 

24 03 37.49 

-15 

C-5 

I 

POTDAM* 

9429 

52 22 55.00 

13 04 01.00 

111 

European 




52 22 52.33 

13 03 55.80 

106 

C-5 

I 

ZVENIG* 

9430 

55 41 37.70 

36 46 03.00 

145 

European 

1 



55 41 36.17 

36 46 00.17 

114 

C-5 

j 


These SAO station positions were derived by using the weighting scheme described in Reference 
2, Section 2. 
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Table in 

STADAN — Optical — Source B 


1 

Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

1BPOIN 

1021 

38°25 '49763 

282°54 '48723 

5 

N.A. 



38 25 49.44 

282 54 48.65 

-50 

C-5 

1FTMYR 

1022 

26 32 51.89 

278 08 03.93 

19 

N.A. 



26 32 53.08 

278 08 03.80 

-42 

C-5 

lOOMER 

1024 

-31 23 30.07 

136 52 11.05 

152 

Australian 



-31 23 26.96 

136 52 14.25 

148 

C-5 

1QUITO 

1025 

-0 37 28.00 

281 25 14.81 

3649 

N.A. 



-0 37 22.63 

281 25 15.23 

3554 

C-5 

1LIMAP 

1026 

-11 46 44.43 

282 50 58.23 

155 

N.A. 



-11 46 37.56 

282 50 58.86 

34 

C-5 

1SATAG 

1028 

-33 09 07.66 

289 19 51.35 

922 

N.A. 



-33 08 58.76 

289 19 52.59 

705 

C-5 

1MOJAV 

1030 

35 19 48.09 

243 06 02.73 

905 

N.A. 



35 19 47.57 

243 05 59.18 

874 

C-5 

1JOBUR 

1031 

-25 52 58.86 

27 42 27.93 

1530 

ARC (Cape) 



-25 53 02.70 

27 42 25.41 

1546 

C-5 

1NEWFL 

1032 

47 44 29.74 

307 16 43.37 

104 

N.A. 



47 44 28.73 

307 16 46.67 

58 

C-5 

1COLEG 

1033 

64 52 19.72 

212 09 47.17 

162 

N.A. 



64 52 17.78 

212 09 37.29 

139 

C-5 

1GFORK 

1034 

48 01 21.40 

262 59 21.56 

253 

N.A. 



48 01 20.81 

262 59 19.55 

200 

C-5 

1WNKFL 

1035 

51 26 44.12 

359 18 14.62 

62 

European 



51 26 40.67 

359 18 08.35 

76 

C-5 


C-5 



Table III (Continued) 


Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

1ROSMA 

1042 

35°12 , 0®.'93 

277°07'41'.'01 

914 

N.A. 



35 12 07.03 

277 07 40.81 

857 

C-5 

1TANAN 

1043 

-19 00 27.09 

47 18 00.46 

1377 

Tananarive 



-19 00 33.26 

47 17 58.89 

1355 

C-5 


Table IV 

STADAN - R/R - Source B 


Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

CAR VON 

1152 

-24°54'14'. , 85 

113°42 , 5S. , 05 

38 

Australian 



-24 54 12.29 

113 42 58.54 

10 

C-5 

ROSRAN 

1126 

35 11 45.05 

277 07 26.23 

880 

N.A. 



35 11 45.15 

277 07 26.02 

823 

C-5 

MADGAR 

1122 

-19 01 13.32 

47 18 09.45 

1403 

Tananarive 



-19 01 19.41 

47 18 07.96 

1382 

C-5 


C-6 




Table V 


NAVY TRANET — Doppler — Source C 


Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

LASHAM 

2006 

51°11 , l(y.'62 

358°58'30.'51 

182 

European 



51 11 07.12 

358 58 24.25 

196 

C-5 

SANHES 

2008 

-23 13 01.74 

314 07 50.59 

608* 

Correga 

Alegre 



-23 13 01,74 

314 07 50.59 

608 

C-5 

PHILIP 

2011 

14 58 57.79 

120 04 25.98 

8 

Tokyo 



14 59 16.42 

120 04 21.61 

-70 

C-5 

SMTHFD 

2012 

-34 40 31.31 

138 39 12.39 

39 

Australian 



-34 40 28.16 

138 39 15.66 

31 

C-5 

MISAWA 

2013 

40 43 04.63 

141 20 04.69 

-10 

Tokyo 



40 43 14.63 

141 19 51.45 

38 

C-5 

ANCHOR 

2014 

61 17 01.98 

210 10 37.46 

61 

N.A. 



61 16 59.60 

210 10 28.60 

44 

C-5 

TAFUNA 

2017 

-14 19 50.19 

189 17 13.96 

6* 

USGS 

1962 Astro 



-14 19 50.19 

189 17 13.96 

6 

C-5 

THULEG 

2018 

76 32 18.62 

291 13 46.72 

43 

N.A. 



76 32 20.72 

291 13 51.07 

-7 

C-5 

MCMRDO 

2019 

-77 50 51.00 

166 40 25.00 

-43 

Mercury 



-77 50 50.58 

166 40 35.02 

-29 

C-5 

WAHIWA 

2100 

21 31 26.86 

202 00 00.63 

380 

Old Hawaiian 



21 31 14.95 

202 00 09.83 

368 

C-5 

LACRES 

2103 

32 16 43.75 

253 14 48.25 

1201 

N.A. 



32 16 43.91 

253 14 45.34 

1162 

C-5 


*MSL 


C-7 





Table \ 

7 (Continued) 



Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

LASHM2 

2106 

51°11 T2V32 

358°58'30 , .'21 

187 

European 



51 11 08.82 

358 58 23.95 

201 

C-5 

APLMND 

2111 

39 09 47.83 

283 06 11.07 

146 

N.A. 



39 09 47.60 

283 06 11.52 

90 

C-5 

PRETOR 

2115 

-25 56 46.09 

28 20 53.00 

1417 

European 



-25 56 49.97 

28 20 50.67 

1595 

C-5 

SHEMYA 

2739 

52 43 01.52 

174 06 51.43 

44 

N.A. 



52 42 56.52 

174 06 44.17 

89 

C-5 

BELTSV 

2742 

39 01 39.46 

283 10 27.25 

50 

N.A. 



39 01 39.23 

283 10 27.72 

-5 

C-5 

STNVIL 

2745 

i 

33 25 31.57 

269 09 10.70 

44 

N.A. 



33 25 31.76 

269 09 09.66 

-10 

C-5 




Table VI 

AIR FORCE — Optical — Source I 


Source 

Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 


ANTIGA 

3106 

17°08 , 51'.'68 

298°12'3T.’41 

7 

N.A. 




17 08 53.88 

298 12 39.19 

-42 

C-5 

E 

GRNVLE 

3333 

33 28 48.97 

268 59 49.17 

45 

N.A. 




33 28 49.15 

268 59 48.12 

-9 

C-5 


GRVILL 

3334 

33 25 31.95 

269 05 11.35 

43 

N.A. 




33 25 32.14 

269 05 10.30 

-10 



USAFAC 

3400 

39 00 22.44 

255 07 01.01 

2191 





39 00 21.99 

255 06 58.32 

2147 


E 

BEDFRD 

3401 

42 27 17.53 

288 43 35.03 

88 

N.A. 




42 27 17.06 

288 43 36.14 

33 

C-5 

E 

SEMMES 

3402 

30 46 49.35 

271 44 52.37 

79 

N.A. 




30 46 49.85 

271 44 51.64 

23 

C-5 


SWANIS 

3404 

17 24 16.57 

276 03 29.87 

83 

N.A. 




17 24 18.90 

276 03 29.71 

18 

C-5 


GRDTRK 

3405 

21 25 47.05 

288 51 14.03 

7 

N.A. 




21 25 48.69 

288 51 15.03 

-48 

C-5 


CURACO 

3406 

12 05 22.11 

291 09 43.76 

23 

N.A. 




12 05 25.49 

291 09 45.16 

-34 

C-5 


TRNDAD 

3407 

10 44 32.78 

298 23 23.67 

269 

N.A. 




10 44 36.16 

298 23 25.43 

210 

C-5 


GRANFK 

3451 

47 56 38.63 

262 37 11.21 

296 

N.A. 




47 56 38.03 

262 37 09.15 

242 

C-5 


TWINOK 

3452 

36 07 25.69 

262 47 04.48 

312 

N.A. 




36 07 25.58 

262 47 02.68 

262 

C-5 


C-9 













Table VI (Continued) 


Source Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

ROTHGR 

3453 

51°25 , 00'.'00 

9 ° 30' 06100 

351 

European 



51 24 57.05 

9 30 00.58 

352 

C-5 

ATHNGR 

3463 

37 53 30.00 

23 44 30.00 

16 

European 



37 53 26.07 

23 44 26.73 

23 

C-5 

TORRSP 

3464 

40 29 18.53 

356 34 41.24 

588 

European 



40 29 14.10 

356 34 36.06 

635 

C-5 

CHOFUJ 

3465 

35 39 57.00 

139 32 12.00 

49 

Tokyo 



35 40 08.96 

139 32 00.19 

75 

C-5 

KINDLY 

3471 

32 22 57.30 

295 19 00.46 

26 

N.A. 



32 22 57.41 

295 19 02.09 

-23 

C-5 

E HUNTER 

3648 

32 00 05.87 

278 50 46.36 

17 

N.A. 



32 00 06.32 

278 50 46.32 

-40 

C-5 

JUPRAF 

3649 

27 01 14.80 

279 53 13.72 

26 

N.A. 



27 01 16.02 

279 53 13.72 

-37 

C-5 

E ABERDN 

3657 

39 28 18.97 

283 55 44.56 

4 

N.A. 



39 28 18.71 

283 55 45.10 

-51 

C-5 

E HOMEST 

3861 

25 30 24.69 

279 36 42.69 

18 

N.A. 



25 30 26.02 

279 36 42.70 

-44 

C-5 

CHYWYN 

3902 

41 07 59.20 

255 08 02.65 

1890 

N.A. 



41 07 58.61 

255 07 59.94 

1845 

C-5 






Table VII 


ARMY MAP SERVICE - SECOR - Source H 



Longitude 


Geodetic 

Height 

(meters) 


Datum 


HERNDN 

5001 

38°59'37:69 

282 o 40'K?.'68 

119 

N.A. 



38 59 37.47 

282 40 17.08 

64 

C-5 

CUBCAL 

5200 

32 48 00.00 

242 52 00.00 

101 

N.A. 



32 47 59.74 

242 51 56.55 

71 

C-5 

LARSON 

5201 

47 11 00.00 

240 40 00.00 

354 

N.A. 



47 10 58.76 

240 39 55.68 

319 

C-5 

WRGTON 

5202 

43 39 00.00 

264 25 00.00 

481 

N.A. 



43 38 59.49 

264 24 58.27 

428 

C-5 

GREENV 

5333 

33 25 32.34 

269 05 10.78 

43 

N.A. 



33 25 32.53 

269 05 09.73 

-10 

C-5 

TRUKIS 

5401 

7 27 39.30 

151 50 31.28 

5* 

Navy Iben 
Astro 1947 



7 27 39.30 

151 50 31.28 

5 

C-5 

SWALLO 

5402 

-10 18 21.42 

166 17 56.79 

9* 

1966 SECOR 
Astro 



-10 18 21.42 

166 17 56.79 

9 

C-5 

KUSAIE 

5403 

5 17 44.43 

163 01 29.88 

7* 

Astro 1962, 
65, Allen 
Sodano Lt 



5 17 44.43 

163 01 29.88 

7 

C-5 

GIZZOO 

5404 

-8 05 40.58 

156 49 24.82 

49* 

Provisional 

DOS 



-8 05 40.58 

156 49 24.82 

49 

C-5 









Table VII (Continued) 


Source 

Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 


TARAWA 

5405 

l 0 21'42’:i3 

172°55'4T.'26 

7* 

1966 SECOR 
Astro 




1 21 42.13 

172 55 47 26 

7 

C-5 


NANDIS 

5406 

-17 45 31.01 

177 27 02.83 

17* 

Viti 

Levu 1916 




-17 45 31.01 

177 27 02.83 

17 

C-5 


CANTON 

5407 

-2 46 28.90 

188 16 43.47 

6* 

1966 Canton 
Astro 




-2 46 28.90 

188 16 43.47 

6 

C-5 


JONSTN 

5408 

16 43 51.68 

190 28 41.55 

6* 

Johnston 
Island 1961 




16 43 51.68 

190 28 41.55 

6 

C-5 


MIDWAY 

5410 

28 12 32.06 

182 37 49.53 

6 

Midway 
Astro 1961 




28 12 32.06 

182 37 49.53 

6 

C-5 


MAUIHI 

5411 

20 49 37.00 

203 31 52.77 

32 

Old Hawaiian 




20 49 25.14 

203 32 01.88 

31 

C-5 

G 

FT WART 

5648 

31 55 18.41 

278 26 00.26 

29 

N.A. 




31 55 18.86 

278 26 00.18 

-27 

C-5 

G 

HNTAFB 

5649 

32 00 04.04 

278 50 43.17 

27 

N.A. 




32 00 04.49 

278 50 43.13 

-30 

C-5 

G 

HOMEFL 

5861 

25 29 21.18 

279 37 39.35 

18 

N.A. 




25 29 22.51 

279 37 39.37 

-44 

C-5 


*MSL 


C-12 



Table VIE 

USC&GS — Optical — Source F 


Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

BELTVL 

6002 

39°01 f 39':03 

283°10 , 26 , . , 94 

45 

N.A. 



39 01 38.80 

283 10 27.40 

-10 

C-5 

ASTRMD 

6100 

39 01 39.72 

283 10 27.83 

45 

N.A. 



39 01 39.49 

283 10 28.29 

-10 

C-5 

TIMINS 

6113 

48 33 56.17 

278 37 44.54 

290 

N.A. 



48 33 55.70 

278 37 44.49 

232 

C-5 


Table IX 

SPEOPT — Optical — Source B 


Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

1UNDAK 

7034 

48°01’21'.'40 

262°59'2T.'56 

255 

N.A. 



48 01 20.81 

262 59 19.55 

201 

C-5 

1EDINB 

7036 

26 22 45.44 

261 40 09.03 

67 

N.A. 



26 22 46.35 

261 40 07.34 

15 

C-5 

1COLBA 

7037 

38 53 36.07 

267 47 42.12 

271 

N.A. 



38 53 35.81 

267 47 40.85 

218 

C-5 

1BERMD 

7039 

32 21 48.83 

295 20 32.56 

21 

N.A. 



32 21 48.94 

295 20 34.18 

-28 

C-5 

1 PUR IO 

7040 

18 15 26.22 

294 00 22.17 

58 

N.A. 


i 

18 15 28.30 

294 00 23.63 

5 

C-5 


C-13 





Table IX (Continued) 


Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

1GSFCP 

7043 

39°01'15'.'01 

283°10'19.'93 

54 

N.A. 



39 01 14.78 

283 10 20.39 

-1 

C-5 

1CKVLE 

7044 

38 22 12.50 

274 21 16.81 

187 

N.A. 



38 22 12.33 

274 21 16.28 

131 

C-5 

1DENVR 

7045 

39 38 48.03 

255 23 41.19 

1796 

N.A. 



39 38 47.54 

255 23 38.52 

1751 

C-5 

1JUM24 

7071 

27 01 12.77 

279 53 12.31 

25 

N.A. 



27 01 14.00 

279 53 12.30 

-38 

C-5 

1JUM40 

7072 

27 01 13.17 

279 53 12.49 

25 

N.A. 



27 01 14.39 

279 53 12.49 

-38 

C-5 

1JUPC1 

7073 

27 01 13.11 

279 53 12.72 

22 

N.A. 



27 01 14.33 

279 53 12.72 

-41 

C-5 

1JUBC4 

7074 

27 01 13.33 

279 53 12.76 

25 

N.A. 



27 01 14.55 

279 53 12.76 

-38 

C-5 

1SUDBR 

7075 

46 27 20.99 

279 03 10.35 

281 

N.A. 



46 27 20.52 

279 03 10.35 

224 

C-5 

1JAMAC 

7076 

18 04 31.98 

283 11 26.52 

485 

N.A. 



18 04 34.20 

283 11 27.03 

423 

C-5 
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Table X 

SPEOPT — Laser — Source B 


Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 

ROSLAS 

7051 

35°11 '40J6O 

277°07'2(?.'23 

879 

N.A. 



35 11 46.70 

277 07 26.02 

822 

C-5 

GODLAS 

7050 

39 01 13.68 

283 10 18.05 

55 

N.A. 



39 01 13.45 

283 10 18.51 

0 

C-5 


Table XI 

INTERNATIONAL — Optical — Source I 


Source 

Name 

Station 

No. 

Latitude 

Longitude 

Geodetic 

Height 

(meters) 

Datum 


DELFTH 

8009 

52 °00' 09^24 

4°22 , 21'.'23 

23 

European 




52 00 06.12 

4 22 15.30 

28 

C-5 


MALVRN 

8011 

52 08 39.12 

358 01 59.49 

111 

European 




52 08 35.68 

358 01 53.03 

125 

C-5 

D 

ZIMWLD 

8010 

46 52 41.77 

7 27 57.56 

898 

BERNE 




46 52 36.73 

7 27 52.54 

907 

C-5 


C-15 



Table XII 
SAO — optical 


Name 

Station No 

10RGAN 

9001 

10LFAN 

9002 

100MER 

9003 

1 SPAIN 

9004 

1TOKYO 

9005 

1NATOL 

9006 

1QUIPA 

9007 

1SHRAZ 

9008 

1CURAC 

9009 

1JUPTR 

9010 

1VILDO 

9011 

1MAUIO 

9012 

OSLONR 

9426 

AUSBAK 

9023 

NATALB 

9029 

AGASSI 

9050 

COLDLK 

9424 

EDWAFB 

9425 

RIGLAT 

9428 

POTDAM 

9429 

ZVENIG 

9430 


Location 

Organ Pass, New Mexico 
Olifantsfontein, South Africa 
Woomera, Australia 
San Fernando, Spain 
Tokyo, Japan 
Naini Tal, India 
Arequipa, Peru 
Shiraz, Iran 

Curacao, Lesser Antilles 
Jupiter, Florida 
Villa Dolores, Argentina 
Maui, Hawaii 
Oslo, Norway 
Woomera, Australia 
Natal, Brazil 

Cambridge, Massachusetts 
Cold Lake, Alberta 
Edwards AFB, California 
Riga, Latvia 
Potsdam, Germany 
Zvenigorod, Russia 


C-16 



Table XIII 


Name 

1BPOIN 

1FTMYR 

lOOMER 

1QUITO 

1LIMAP 

1 SAT AG 

1MOJAV 

1JOBUR 

1NEWFL 

1COLEG 

1GFORK 

1WNKFL 

1ROSMA 

1TANAN 


Name 

CAR VON 
ROSMAN 
MADGAR 


STADAN - Optical 


Station No. 

Location 

1021 

Blossom Point, Maryland 

1022 

Fort Myers, Florida 

1024 

Woomera, Australia 

1025 

Quito, Ecuador 

1026 

Lima, Peru 

1028 

Santiago, Chile 

1030 

Mojave, California 

1031 

Johannesburg, Union of South Africa 

1032 

St. John's, Newfoundland 

1033 

College, Alaska 

1034 

East Grand Forks, Minnesota 

1035 

Winkfield, England 

1042 

Rosman, North Carolina 

1043 

Tananarive, Madagascar 


Table XIV 
STADAN - R/R 


Station No. 

Location 

1152 

Carnarvon, Australia 

1126 

Rosman, North Carolina 

1122 

Tananarive, Madagascar 
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Table XV 


NAVY TRANE T - Doppler 


Name 

Station No. 

Location 

LASHAM 

2006 

Lasham, England 

SANHES 

2008 

Sao Jose dos Campos, Brazil 

PHILIP 

2011 

San Miquel, Philippines 

SMTHFD 

2012 

Smithfield, Australia 

MISAWA 

2013 

Misawa, Japan 

ANCHOR 

2014 

Anchorage, Alaska 

TAFUNA 

2017 

Tafuna, American Samoa 

THULEG 

2018 

Thule, Greenland 

MCMRDO 

2019 

McMurdo Sound, Antarctica 

WAHIWA 

2100 

South Point, Hawaii 

LAC RES 

2103 

Las Cruces, New Mexico 

LASHM2 

2106 

Lasham, England 

APLMND 

2111 

APL Howard County, Maryland 

PRETOR 

2115 

Pretoria, Union of South Africa 

SHEMYA 

2739 

Shemya Island, Alaska 

BELTSV 

2742 

Beltsville, Maryland 

STNVIL 

2745 

Stoneville, Mississippi 
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Table XVI 


Name 

AN TIG A 

GRNVLE 

GRVILL 

USAFAC 

BEDFRD 

SEMMES 

SWANIS 

GRDTRK 

CURACO 

TRNDAD 

GRANFK 

TWINOK 

ROTHGR 

ATHNGR 

TORRSP 

CHOFUJ 

KINDLY 

HUNTER 

JUPRAF 

ABERDN 

HOMEST 

CHYWYN 


AIR FORCE - Optical 


Station No. 

Location 

3106 

Antigua Island, Lesser Antilles 

3333 

Stoneville, Mississippi 

3334 

Stoneville, Mississippi 

3400 

Colorado Springs, Colorado 

3401 

L. G. Hanscom Field, Massachusetts 

3402 

Semm es Island , Georgia 

3404 

Swan Island, Caribbean Sea 

3405 

Grand Turk, Caicos Islands 

3406 

Curacao, Lesser Antilles 

3407 

Trinidad Island 

3451 

Grand Forks , North Dakota 

3452 

Twin Oaks, Oklahoma 

3453 

Rothwesten, West Germany 

3463 

Athens, Greece 

3464 

Torrejon de Ardoz , Spain 

3465 

Chofu, Japan 

3471 

Kindly AFB, Bermuda 

3648 

Hunter AFB, Georgia 

3649 

Jupiter, Florida 

3657 

Aberdeen, Maryland 

3861 

Homestead AFB, Florida 

3902 

Cheyenne , Wyoming 
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Table XVII 

ARMY MAP SERVICE - Secor 


Name 

Station No. 

Location 

HERNDN 

5001 

Herndon, Virginia 

CUBCAL 

5200 

San Diego, California 

LARSON 

5201 

Moses Lake, Washington 

WRGTON 

5202 

Worthington, Minnesota 

GREENV 

5333 

Greenville, Mississippi 

TRUKIS 

5401 

Truk Island, Caroline Islands 

SWALLO 

5402 

Swallow Island, Santa Cruz Islands 

KUSAIE 

5403 

Kusaie Islands, Caroline Island 

GIZZOO 

5404 

Gizzoo, Gonzongo, Solomon Islands 

TARAWA 

5405 

Tarawa, Gilbert Islands 

NANDIS 

5406 

Nandi, Vitilevu, Fiji Islands 

CANTON 

5407 

Canton Island, Phoenix Islands 

JONSTN 

5408 

Johnston Island, Pacific Ocean 

MIDWAY 

5410 

Eastern Island, Midway Islands 

MAUIHI 

5411 

Maui, Hawaii 

FTWART 

5648 

Fort Stewart, Georgia 

HNTAFB 

5649 

Hunter AFB, Georgia 

HOMEFL 

5861 

Homestead AFB, Florida 
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Table XVIII 
USC & GS — Optical 


Name 

Station No. 

Location 

BELTVL 

6002 

Beltsville, Maryland 

ASTRMD 

6100 

Beltsville, Maryland 

TIM INS 

6113 

Timmins, Ontario 


Name 

Table XIX 
SPEOPT - Optical 

Station No. Location 

1UNDAK 

7034 

Univ. North Dakota, Grand Forks 

1EDINB 

7036 

North Dakota 
Edinburg, Texas 

1COLBA 

7037 

Columbia, Missouri 

1BERMD 

7039 

Bermuda Island 

1PURIO 

7040 

San Juan, Puerto Rico 

1GSFCP 

7043 

GSFC, Greenbelt, Maryland 

1CKVLE 

7044 

Clarksville, Indiana 

1DENVR 

7045 

Denver, Colorado 

1JUM24 

7071 

Jupiter, Florida 

1JUM40 

7072 

Jupiter, Florida 

1JUPC1 

7073 

Jupiter, Florida 

1JUBC4 

7074 

Jupiter, Florida 

1SUDBR 

7075 

Sudbury, Ontario 

1JAMAC 

7076 

Jamaica, B. W. I. 
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Table XX 


SPEOPT — Laser 


Name 

Station No. 

Location 

ROSLAS 

7051 

Rosman, North Carolina 

GOD LAS 

7050 

GSFC, Greenbelt, Maryland 


Table XXI 

INTERNATIONAL - Optical 


Name 

Station No. 

Location 

DELFTH 

8009 

Delft, Holland 

MALVRN 

8011 

Malvern, England 

ZIMWLD 

8010 

Berne, Switzerland 


C-22 


2.0 SOURCES 


The following sources were used to obtain the original datum positions: 

Symbol Source 

A Geodetic Parameters for a Standard Earth Obtained from 

Baker-Nunn Observations; September 1966; Smithsonian 
Astrophysical Observatory. 

B Goddard Directory of Tracking Station Locations; August 

1966; Goddard Space Flight Center. 

C NWL-8 Geodetic Parameters Based on Doppler Satellite 

Observations; July 1967; R. Anderle and S. Smith, Naval 
Weapons Laboratory. 

Since the above official documents did not contain all those positions that 
were to be transformed, it was necessary to contact other sources for the posi- 
tions of the remaining stations. These sources are: 


Symbol 

D 

E 

F 

G 

H 


Source 

Private communication with personnel at SAO; K. Haramundanis; 
E. Miller; A. Girnius. 

Private communication with 1381 Geodetic Survey Squadron, 
USAF; S. Tischler. 

Private communication with personnel at USC&GS; B. Stevens. 

Private communication with personnel at U.S. Army Engineers 
Topographic Laboratories; L. Gambino. 

Private communication with NASA Space Science Data Center; 

J. Johns; D. Tidwell. 

General Station Data Sheet — GEOS-A-Project Manager NASA 
Headquarters. 
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